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Abstract. We are interested in the motion of a classical charge coupled to the 
Maxwell self- field and subject to a uniform external magnetic field, Bq. This is a 
physically relevant, but difficult dynamical problem, to which contributions range 
over more than one hundred years. Specifically, we will study the Sommerfeld-Page 
approximation which assumes an extended charge distribution at small velocities. 
The memory equation is then linear and many details become available. We discuss 
how the friction equation arises in the limit of "small" Bq and contrast this result 
with the standard Taylor expansion resulting in a second order equation for the 
velocity of the charge. 



1 Introduction 



Radiative damping must be taken into account in the design of synchrotron radia- 
tion sources pj. Amongst other tasks, ultra strong lasers are expected to provide a 
quantitative test of radiative damping forces in the ultra relativistic regime P] . For 
stellar jets and other highly accelerated beams of charged particles radiative friction 
plays a dominant role [1]. The ground work on radiative friction was accomplished 
already Abraham [3] and Lorentz |1] , who undertook a vast effort to derive effective 
equations of motion for the charged particle starting from a model of an extended 
charged body coupled to the Maxwell field. Later on fully relativistic models fol- 
lowed. I refer to recent studies [5], E] and to the monographs [3 El E] , in which prior 
work is put in perspective. To be clear, we will consider classical charges coupled 
to the classical electromagnetic field. When the interest is in the radiative decay of, 
say, an excited hydrogen atom, then a quantum version of the theory is required, 
which is beyond the scope of this note. 

Granted a few exceptions, in all derivations of the effective equations of motion 
for the charged particle one uses Taylor expansion in a small parameter, which leads 
to a differential equation of second order in the velocity. In such a procedure it 
is never detailed how to match the true initial data with the initial conditions for 
the effective dynamics. In addition the second order differential equation generates 
spurious solutions which are a mere artifact of the approximation method. The 
case study presented here is so simple that both defects can be easily exhibited 
and remedied. While we discuss the extension to more physical models, we believe 
that a satisfactory derivation of the equations of radiative friction are still to be 
accomplished. 

For the sake of discussion I prefer to fix a very concrete physical set-up, namely 
the motion of a charged particle subject to a uniform external magnetic field and 
coupled to its self-generated Maxwell field. The coupling to the self-field makes our 
dynamical problem difficult. To keep the presentation as simple as possible and 
to be able to provide a complete argument, we will study an oversimplified model 
known as Sommerfeld-Page equation. In the Sommerfeld-Page equation the charge 
distribution of the particle is taken to be uniform over a sphere of radius R and 
normalized to e. The charge distribution is supposed to be rigid and independent 
of velocity in the given inertial frame. Also small velocities are assumed. At the 
initial time t = one has to specify position and velocity of the center of the charge 
distribution and the self-field. One natural choice for the self-field is the Coulomb 
field generated by the charge distribution at rest. An even more physical choice for 
the self-field would be the field as if the particle has been traveling with constant 
velocity Vq for alH < 0. There are other choices and the initial Maxwell self- field is 
to some extent arbitrary. For sure, finite energy and some smoothness is required. 
On the other hand, only as one example, an incident plane wave would have to be 
taken into account explicitly in the Sommerfeld-Page equation. 
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With these assumptions the equations for the motion of the particle read 



Here mo is the mechanical (bare) mass of the particle, e the charge, and c the velocity 
of light. The magnetic field points along the 3-direction with magnitude Bq. The 
velocity is decomposed in the component v orthogonal and the component f 3 parallel 



to the magnetic field. If = (f 1,^72), then v = {—V2,Vi). In the literature (1.1), 



(1.2) is referred to as Sommerfeld-Page equation [ini [H], see also Appendix A for 



more details on the approximations involved. We will discuss only (1.1), since (1.2) 



is automatically included as the special case Bq = in (1.1). 
To streamline the notation, we set 

eBo 2R e2 

u= , T = — , a = - —. (1.3) 

nioc c SirniQC'^T'^ 

Thus u measures the strength of the external magnetic field and a the coupling 



strength to the self-field. (1.1) becomes then 

d 



v{t) = ujv{t)^ + a{v{t - r) - v{t)) . (1.4) 

(JjL 



The natural dimensionless parameters will be ut and ar. (1.4) is a delay equation 
because the current change in velocity depends on the past history. Also the notion 
differential delay equation is very common. In statistical mechanics one uses the 



term memory equation, indicating again the dependence on the past. Eq. (1.4) has 
to be supplemented with the initial condition 

v{t) = u{t), 0<t<T, (1.5) 

with some specified function u{t) encoding the interaction with the self- field up to 
time r, an information which usually is not available in explicit form. In fact we 
will see that the precise details of u{t), except for the average and the end point 
u{t), play no role. 



To obtain from (1.1) an effective differential equation seems to be an easy task. 



Since R is small, one Taylor expands in R up to second order and arrives at 

("^0 + ^) tM^) = — + ^^^W ' (1-6) 



where we use w{t) to distinguish from the true solution v{t). Formally one expects 
that 

v{t) - w{t) = 0{R) . (1.7) 
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The prefactor of w{t) is regarded as effective mass, 



rries = mo + 



(1 + aT)mQ , 



If one really wants i? — )■ 0, following Dirac [12], one takes the double limit i? — )■ 0, 
mo — 7- — oo with rrics > fixed. Then it is claimed that 



lim v(t) = wit) . 

R-s>0 



:i.9) 



To fix a unique solution to (1.6) one has to prescribe ti>(0) and w{fi). On the 
other hand for (1.4) we need to specify a whole initial function w(t), < t < r. A 



complete derivation of (1.6), or (1.9), thus must give a prescription how to link u{t) 



to wifi) and w(Q). This part is missing and the Taylor expansion does not tell us. 
In fact, the situation is even more complicated. The phase space for ( |1.6 ) is 



X . 



^ Embedded is the two-dimensional plane Cuj with G C^. If (i(^(0), w{Q)) G 
Ct^;, then the motion is stable and spirals towards 0. On the other hand for initial 
data off Ctj, the motion diverges exponentially fast to infinity. Since such run away 
solutions have never been observed, I regard them as unphysical. But then one has 
to explain why for a broad class of initial u{t), < t < r, one always has initial 
data for (1.9) which lie in C^^. Again the Taylor expansion does not tell us. 



For the Sommerfeld-Page equation a straightforward alternative to the Taylor 



expansion (1.6) is available. Eq. (1.4) is the standard and most elementary example 



of a linear differential-difference equation which is exhaustively covered in the liter- 
ature [131 [H]. I will use here the conventional method for construction of solutions. 
Given the linear structure, it is rather surprising that so far the Sommerfeld-Page 
equation has been analyzed only rudimentarily [151 HSl El [El [El 120] . 

Our task is to study the solution to (1.4), (1.5). Roughly one expects that 



the synchrotron motion is damped and the particle will come to rest in the long 
time limit. In principle one could explore properties of (1.4) for all u and a. But 
it is of advantage to keep the physical parameter range in mind. The classical 
electron radius is defined by the condition that the mechanical rest energy equals 
the electrostatic energy of the charge distribution. This yields R ~ 10~^^ cm. Then 
ar = 10~^ and for a magnetic field of 10® Tesla one arrives at ut = 10~^. The 
strongest laboratory fields are of order 10 Tesla. Thus it will be safe to eventually 
consider only the regime < ar ^ 1, \ujt\ ^ 1. 

Before entering into a more detailed analysis, let me summarize already now the 
main findings. 

(1) Zero external magnetic field, Bq = 0. 

(i) If rries/mQ > 0, then the motion is stable. On the other hand if rries/mQ < 0, 
generically the solution increases exponentially. Expressed in dependence on mo, if 
mo > 0, then nies > 0, and the motion is stable. If — e^/37rc^r^ < mo < 0, then the 
motion is unstable but still nies > 0. If mo < —e"^ /Sttc^t^, then nies < and the 
motion is again stable. 
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(a) If mofj/mo > 0, then v{t) v exponentially fast as t — )■ oo, where 

V = {1 + aTy\u{T) + a [ dsu(s)) , (1.10) 

Jo 

to say, there is a decay constant '-far < 1 such that 

\vinr)-v\<Coi7arr (1-11) 

with integer n > and some constant cq. At ar = 0.278, the decay constant 
7a,- ~ 0.07 and decreases rapidly as ar decreases to 0. At ar = —1, the decay 
constant •jaT ~ 0.25. 

(2) Non-zero external magnetic field, Bq ^ 0, \ujt\ ^ 1, mg > 0. 

(Hi) If mo > 0, |ci;r| ^ 1, and ur ^ 0, then v{t) — )■ exponentially fast as t — )■ cxo 
with an upper bound 

\v{nr)\<cJl-^-^^i-^fy . (1.12) 



(iv) For mo > and at fixed R, if \u)t\ ^ 1, the solution to (1.4) is well approximated 
by 

(1 + aT)^w(t) = Lowit)^ - luJ^Tarn + ary^wit) . (1.13) 
at ^ 

More precisely, one can choose an initial condition w{0), such that w{0) = v+0{ut) 
and 

\v{t)-w{t)\<0{{tuTf) (1.14) 

for all t/r > | log(t<;r)^/ log7Q,,-|. Note that there is a slip in the initial velocity 
to V + 0{ut). In the physical regime, instead of the more complicated memory 



units of ( 1.1 ) reads 



equation (1.1), one can work with the friction differential equation (1.13), which in 



d eBo 

meE—v[t) = 

dt c 



, , I / eBo \2 , , 



:i,i5) 



A detailed discussion and derivation of these results is presented in Sections 2 
and 3, where we also explain the point charge limit i? — >■ 0. In our context Landau 



and Lifshitz [21] provide a simple recipe how to proceed from (1.6) to arrive directly 



at (1.15). This point will be discussed in more detail in Section 4. Based on the 



novel insights from the Sommerfeld-Page equation we return to the full model, a 
charge coupled to its Maxwell self-field, and summarize of what has been achieved 
and what still has to be done, in our opinion. 
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2 Stability of the memory equation 



For the study of the Sommerfeld-Page equation it will be convenient to use complex 
notation. Then v{t) = Vi{t) + iv2{t) G C and u{t) = Ui{t) + iu2{t) G C. With this 
notation (1.4) reads 



d 
dt 



and, according to (1.5), 



;(t) = -iujv{t) + a{v{t - r) - v{t)) 



v{t) = u{t) for < t < r. 



To find out v{t + t), we integrate (2.1) on both sides. Then 

v{t + t)= e-("+^'^)*M(r) + a f dse-'-''+'''^^'-'^u{s) = {K^u){t) 



(2.1) 



(2.2) 



(2.3) 



< t < r. We regard (2.3) as a linear map which out of v in the interval [0, r] 
makes v in the interval [r, 2r]. We will call K^^ the Sommerfeld-Page operator. Note 
that is a non-symmetric operator. More details are explained in Appendix B. 
To obtain the solution v{t) for all times one has to simply iterate the map (2.3) with 
the result 

v{nT + t) = {{K^ru){t), (2.4) 

n integer and < t < r. Clearly, {Ki^u){0) = u{t). If u is arbitrary, then v{t + r), 
< t < r, is continuous and so is v{t) for all r < t < oo by iteration. 

The stability of solutions is determined by the eigenvalue problem for K^, 



KujU = Xu - 



(2.5) 



Stable dynamics is equivalent to |A| < 1 for all eigenvalues. Differentiating (2.5) 
with respect to t yields 



X^uit) = au(t) — X(a + iuj)u(t) , 
dt 

which implies that the eigenvector u is an exponential function. Setting 



one concludes 



u{t) 



X 



zt/r 



and |A| < 1 precisely ii <0. Inserting in (2.6) one obtains 

z = —iuT + ar(e~^ — 1) . 



(2.6) 

(2.7) 
(2.8) 

(2.9) 
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We will study case by case and, to simplify notation, switch to dimensionless vari- 
ables. Since we want to allow mo < 0, it is convenient to set r = 1 for a while. 



(i) The case u = 0. The eigenvalues of Kq are determined by 

-1). 



z = a[e 



(2.10) 



We set z = X + iy, x,?/ G M, and first study solutions for which Qz = 0. They are 
determined by 

x + a = ae-''. (2.11) 

X = = xo{a) is always a solution. Let us introduce the two critical points, a+, 
a- = —1 and a;+ is determined by 



a , e"+ = e"^ 



0.278. 



(2.12) 



(2.11) has one further solution, xi{a), provided a < 0, where xi{a) < for < 
a < and xi{a) > for a < a_. Thus the motion is unstable for 171^^/1710 < 0. 



Equating absolute value and phase of (2.10), one arrives at 



(A) 
(P) 



(x + af + 



X + a = —ycoty . 



(2.13) 
(2.14) 



The intersections of the curves (A) and (P) define the additional eigenvalues. (P) 
has branches labeled by m G Z corresponding to (m — l)7r < y < mvr, m < — 1, 
— TT < y < n, rriTi < y < {m + l)7r, m > 0. The 0-th branch has its apex at 
X = —1 — a, y = and converges to ±7r as x — )■ 00. The other branches are strictly 
increasing for m > 0, strictly decreasing for m < 0, and cover the full real line, see 
Fig.l. 











n 










1 " 















Figure 1: Curve (P) for a = 0. Eigenvalues result by intersection with one of the 
curves from Fig. 2, resp. Fig. 3. 



Next we discuss (A), which has the following solutions at y = 0, see Fig. 2 and 
3: Xo(a) = is a solution for all a. Xi{a) and X2(a) exists only for a < «+. One 
has xi{a) < ^2(0;) < for < a < a+, xi{a) < < ^2(0;) for a_ < a < 0, and 
< Xi{a) < X2{<y) for a < a_. (A) has always a branch which diverges to ±00 
as X —7- —00. Its apex is located at (xi(a;),0) for a_ < a < a;+ and at (0,0) for 
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Figure 2: (a) the curve (A) for a = 0.27835 < a+ and (b) for a = 0.279 > a+. • is 
an eigenvalue. 




a < a_, a+ < a. In addition (A) has a bubble with apices at (x2(a),0) and (0,0) 
for and at {xi{a), 0) and (0:2(0;), 0) for a < a_. 

The 0-th branch of (P) never intersects (A). The branches m 7^ of (P) intersect 
the diverging branch of (A), which thus yields only stable eigenvalues. For < a < 
a+ they satisfy < —1 — a+. The spectrum of eigenvalues of at a; = is 
plotted schematically in Fig. 4. The decay constant •yar is the radius of the smallest 
disk which contains all eigenvalues except for the eigenvalue 1. For a < a_ the 
intersection points lie to the left but very close to the imaginary axis. Hence, the 



stability of motion is determined by the solutions to (2.11). 

(ii) The case w 7^ 0, ^ 1. The eigenvalues are determined by 

z = -icu + a{e'' - 1) , (2.15) 

which induces a small shift. Of interest is the shift of the eigenvalue A = 1, i.e. 
z = 0. Expanding at z = 0, one arrives at 

z = -iu + a{-z + Iz^) , (2.16) 

which for small u has the solution 

z = -^^—u^-i-^—u. (2.17) 

2(1 + aY (1 + a) ^ ' 
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Figure 4: Schematic spectrum of K^. • are eigenvalues for oj — 0. The radius of the 
shaded disk is the decay constant, o indicates the shift of the eigenvalues by taking 

For a > the motion is stable, while for —1 < a < it is unstable. At a = — 1 there 
is a double degenerate eigenvalue 0, which splits into the solutions of z"^ = —2iu). 
Hence the unstable motion persists for a < —1. For <^ 1 and a > 0, there is a 
single maximal eigenvalue, Amax, with the property 

\X^l^l-l^J^r + 0(.'). (2.18) 

which implies a slow inward motion. 

3 Long time behavior and comparison dynamics 

These two topics are not specific for a charge coupled to the Maxwell field. The long 
time behavior is a first central task in the analysis of a dynamical system. Compar- 
ison dynamics is the issue of approximating the full dynamics through simplified, 
more tractable equations and typically requires to have a small parameter at one's 
disposal. The generic form of such a comparison is to first match the initial data 
then to provide error bars which specify how well and over what time scale a solution 
of the true dynamics is shadowed by the comparison dynamics. In the following we 
denote by v{t) the true dynamics, as before, and the approximation by w{t), both 
taking values in C. 

3.1 Long time asymptotics 

To determine the long time behavior we need the maximal eigenvector and the 
spectral gap. As explained in more detail in the appendix, the natural space for the 
Sommerfeld-Page operator is the Hilbert space = C © L^([0, r], dt). A vector 
{uo,u) e H consists of the constant uq = u{t) and the function u{t), < t < t. 
The spectrum of K^^ acting on H is denoted by a{Ki^). As shown in Appendix B, 
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K^^ has the family of right eigenvectors \gx) and of left eigenvectors (/a|, A G a{Ki^). 
such that 

{K^r= A"(1 + A-W)-1|(7a)(/a|. (3.1) 



Explicit formulae for \g\)., \f\) are given in (B.20) together with (B.15), (B.16) at 
q = \. 

We assume now mo > and w = 0. Then 1 is the unique maximal eigenvalue 
with left eigenvector \gi) = (1,1) and right eigenvector = {I, a). Hence, with 
Pi = |l + a)-i|^?i)(/i|, 



{Piu){t) = (1 + ar)"^(M(r) + a I dsu{s)) = v 

Jo 



(3.2) 



for all t, < t < r, which equals Eq. ( 1.10 ) in complex rather than 2- vector notation 



The estimates for the decay constant are derived from the graphical representation 
of Section 2. 

Turning to the case u ^ 0, the maximal eigenvalue Amax has the property that 
|Amax| < 1- Therefore v(t) decays to as t — >■ oo with decay constant jAmaxI- This 
observation yields ( 1.12[ ). 



3.2 Approximation of small 5o-field 

We consider the case \ujt\ ^ 1, which physically corresponds to a "weak" magnetic 
field and assume rriQ > 0, i.e. a > 0. Denoting the integer part by [■], it holds 



(3.3) 



The maximal eigenvalue of is given by the solution to ( |2.15 ) close to z = 0. The 
decay constant is ■jar < 0.07. Therefore, with A = Amax, 

v{t) = A[*/^l(l + X-'aT)-'\g^){h\u){t - [t/r]r) + 0{{j^rf/^^) . (3.4) 

For \ujt\ ^ 1, one obtains 

exp [( - i(l + ar) - a;r(2(l + ar)^) ^w^r)t] + 0((cur)^) . (3.5) 

The initial value is v + 0{ujt). Finally, if t/r > | log(u;r)^/ log7Q,r|, then (7ar)^*^''^ < 
(wr)^. Thereby we have established (1.13) and (1.14). 



3.3 The point charge limit 

Following Dirac [I2] we consider the limit 



r — > for fixed rrics = ttiq + 



37rc'^ 



> 0. 



r 



(3.6) 
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Then, for small r, 



and 



ar = —1 — ar , a = — ^wigg 



UT = —br'^ , b 



c 



(3.7) 



(3.8) 



Since ar < —1, but close to —1, we now have to take into account the two maximal 
eigenvalues, A+, A_, close to 1. The eigenvalue equation reads 



z + iuT = aT{e~^ — 1) . 

Setting z = KT one finds 

A± = 1 + TK± , K± = — a(l ± a/1 + i26/a) . 
In the limit r — )■ it holds 

(l + r/t± )[*/"] ^e^^\ 



(3.9) 



(3.10) 



(3.11) 



In differential form, (3.11) agrees with Lorentz-Dirac equation (1.15). 



Somewhat more tricky is the computation of the initial conditions. In the sum 



(|3.1|) only the eigenvalues X± contribute. All other eigenvalues satisfy |A| < 0.25, 

t < r, 

(3.12) 



independent of r. We set A = 1+tk and omit the ± index. For given u(t), <t < t, 
one has to study 

\gx){fx\u)il + aT\-Y' 



in the limit r — )■ 0. 



Since, according to (B.16), (B.20), 



\gx) = {l,z ^exp[at(2; ^ - 1) - iut]) 



with A = e^, one finds 



lim I^a) = (1, 1) 



The left eigenvector is given by 

l/A) = (l,/A(t)) 

with A = e^ and 



(3.13) 
(3.14) 

(3.15) 



f^(t) = -(1 - 2itu{z/a)) '2ia;e-("~'")("-*) 

+ (l - 2iuj{z/a))~\a/z) exp[at{z-^ - 1) - iw(^ - t)] . (3.16) 



Hence the scalar product becomes 

{fx\u) = n(r) + / dsuis)fxisy 



(3.17) 
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The contribution from the first summand of (3.16) is 0{t^). For the second term we 
expand as u(t) = u{0) +^'(0)^ and obtain the scalar product of (3.17) as ^t{u'{0) — 
K.u{0)). Hence 



\im \gx){fx\u){l + arX''y' = {2{K-a)) \u'{0) - ku{0)) 



r-5-O 



(3.18) 



We conclude that in the point charge limit v{t) w{t) a.s t ^ constrained by 

(3.19) 



(3.6) and the limit motion is given by 



a=± 



is the stable and k_ the unstable mode. Unless there is a specific relation 



between m(0) and 'u'(O) the motion (3.19) is unstable. Dirac postulated his asymp- 
totic condition, which in our context would single out the stable mode. From the 
point of view of the memory equation, as soon as mo < with rngg > 0, the motion 
is unstable. Thus it is not surprising that also the r — )■ limiting motion is un- 
stable. Of course, one could adjust the initial conditions to be orthogonal to |/a_), 
thus maintaining a stable solution throughout the limiting procedure. Clearly, such 
a constraint would have to be justified by some additional reasoning beyond the 
Sommerfeld-Page equation. 



4 The Landau-Lifshitz approximation 

In the Landau-Lifshitz approximation the external -Bo-field is considered to be small. 



To introduce the small parameter explicitly we substitute in (1.4) u by ecu with 
< e ^ 1 a dimensionless parameter. To make up for the small field one has to 
consider times of order e~^. The rescaled time, again denoted by t, is in units of 



single revolutions. With both changes (1.4) becomes 



e^v{t) = euv{t)^ + a{v{t - er) - v{t)) . (4.1) 
Taylor expanding in e yields 

d (f 
(1 + ar)—v{t) = uv{t)^ + \aT'^e—v{t) . (4.2) 

The small parameter e appears now in front of the highest derivative, which is the 
hallmark of singular perturbation theory [22] • The center manifold to order e is 



singled out by re-substituting the leading order in (4.2) with the result 



;i + ar)^v{t) = uv{t)^ - \ea{-^fv{t) ■ (4-3) 
at i + ar 
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It seems that by a clever short-cut we have directly arrived at (1.13). If one 
compares with (1.6), then (4.3) corresponds to the weak i?o-field approximation of 
the Lorentz-Dirac equation restricted to the stable plane C^^ of physical solutions. 

If one tries to control the error in the step from (4.1) to (4.2), one runs into the 
same difficulty as with the Lorentz-Dirac equation. Unstable solutions increase as 
e*/^ and it remains mysterious by which fine tuning of the initial conditions such 
runaways are avoided. But since we know from Section 3 that (4.3) is a controlled 
approximation to the true solution, we conclude that Taylor expanding (4.1) gener- 
ates spurious solutions. They are an artifact of the approximation scheme used and 
not an intrinsic property of the Sommerfeld-Page equation. 



5 Discussions and conclusions 

Physically one would expect that the results from the case study are of a much 
greater generality. Of course the goal would be to treat a fully relativistic matter 
field coupled to the Maxwell field. Not so many rigorous results are available then 
and we propose to first have a look at the semi-relativistic Abraham model. 

The dynamical variables are position q, the velocity v, and the self-fields E{x, t), 
B{x,t). In addition there is a uniform external magnetic field Bq. Then the coupled 
equations of motion read 

c~'-B = -VxE, c-'-E = VxB- -^ip{- - q)v , (5.1) 
V-E = eip{--q), V-B = 0, (5.2) 

—mb^vit) = -v(t) X Bo + - / dx(f(x - q(t)) (E(x, t) + v(t) x B(x, t)) . (5.3) 
at c c J ^ 

Here ip is the rigid charge distribution, (p > 0, ip radial and smooth, (fix) = for 
l^^l ^ J dxip{x) = 1. TTT-b is the mechanical bare mass and 7(t>) = (1 — (i>/c)^)~^/^. 
The initial position and velocity for the particle is ^q, Vq, \vq/c\ < 1. For the self-field 
there is some freedom, as discussed already. To be specific we set 

E{x,0) = E,, {x - go) , B{x, 0) = B^, {x - q^) , (5.4) 

where -E^^, B^^ are the minimizers of the total, i.e. field plus particle, energy at 
given qQ,vo. For an exphcit formula we refer to [9], Eqs. (4.5) to (4.7). 

As for the Sommerfeld-Page equation, the motion in the plane orthogonal to Bq 
can be decoupled. For notational simplicity, let us set Bq = (0,0,i?o)- We assume 

g3(0) = 0, vsiO) = (5.5) 

and require for the initial fields that 

E3{xi,X2, Xs) = -E3{xi,X2, -X3) , Bj{xi, X2, X3) = -Bj{xi, X2, -X3) , j = 1, 2 . 

(5.6) 
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Then from (5.1) - (5.3) it follows that these properties remain valid for all t. Thus 



for initial fields with the particular symmetry (5.6) the motion of the charge stays 
in the 1-2 plane. 

We now follow the scheme of Section 4. In the fully relativistic setting this yields 
the effective equations of motion 

valid for "small" Bq. Here tuq is the rest mass of the particle. If one would aim at 



a precision of order |-Bo|, then the initial conditions for (5.7) are simply Qq, Vq. 

The semi-relativistic character of Abraham model is merely reflected in a velocity 
dependent effective mass which differs from the one of a relativistic particle. The 
relevant computation can be found in [9], Chapter 4.1. With this difference one can 
still work out the re-substitution as indicated in Section 4. The expressions become 
unwieldy and there is no particular reason to display them here. The structure of 



the effective equations of motion is similar to (5.7). 

For our particular initial condition, if Bo = 0, then the particle would simply 
continue to travel with velocity Vq. Thus one might guess that the effective equations 



of motion (5.7) should be solved with initial conditions v{0) = Vq. As discussed at 
length in our case study, to have an approximation which includes friction and thus 
allows for an error 0{BI) only, one expects that the initial conditions must be 
corrected slightly by a term of order |-Bo|- 

The details of such a program are carried out in [23]. As emphasized already, 
upon Taylor expanding one arrives at a comparison dynamics which has runaway 
solutions making difficult the precise control of errors. The best we could accomplish 
is an estimate of the difference in energies of the true and comparison dynamics. A 
pointwise comparison of the trajectories is missing. 

More tractable is the issue of the long time asymptotics. Under certain condi- 
tions it is proved that, in case Bq = 0, limt_>oo v{t) = v [211 123 ES]- Also, if the 
charge is subject to an external electrostatic potential, then in for long times the 
particle comes to rest at a critical point of the potential [27]. On the other side, 
for a uniform external magnetic field Bq ^ the energy balance used in [27] is not 
strong enough to control the long time behavior and it remains as an open problem 
to establish that the charge comes to rest in the limit t — )■ oo. 
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A Appendix: Derivation of the Sommerfeld-Page 
equation 

In standard derivations of the Sommerfeld-Page equation, one expands the solution 
of the full Abraham model in i? as a small parameter retaining only the terms 
linear in v. In jS], Chapter 7, we provide a somewhat different perspective. No 
originality is claimed. Starting from the full Abraham model one can write down 
an exact memory equation for the motion of the charge. There are two terms: one 
term results from the initial conditions of the Maxwell field and a memory term, 
which contains only the self-generated field. The charge distribution is assumed to 
be radially symmetric and supported in a ball of radius R, but otherwise arbitrary, 
and the Maxwell field is assumed to be co-moving. Then memory refers to times 
up to 2R/c into the past and after a time of order 2R/c the initial Maxwell field 
no longer influences the motion of the charge. No expansion is employed. If one 
now specializes to a charge distribution uniformly over a sphere of radius R and 
retains in the memory term only contributions linear in v, then one arrives at the 



Sommerfeld-Page equation (1.1). 



B Appendix: The Sommerfeld-Page operator 

We provide a more detailed analysis of the Sommerfeld-Page operator which 
reads 

{K^u){t) = e-("+''^)*M(r) + a [ dse-^''+''^^^'-'\{s) , < t < r . (B.l) 

Jo 



The eigensolutions to (B.l) 



K^ux = Xux , (B.2) 

are of the form 

U),{t) = e"*/" , A = e^ (B.3) 

Here A G <y{K), the spectrum of K, and A is implicitly determined through the 
solutions of 

z + iojT = ar{e''' - 1) , (B.4) 



which coincides with (2.9). Assuming the initial condition u as the linear superpo- 
sition 

u{t)= caMaW, (B.5) 

one finds 

{K^)Mt)= Yl CAA"MA(t). (B.6) 
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For the particular case w = and 1 + ar > 0, A = 1 is the unique maximal 
eigenvalue. We then split as 



{Ko)Mt) = ci+ CAA"MA(t), (B.7) 

a{K) = a{K) \ {1}. Therefore 

liKoTuit) - ci\ < { Yl IcaDt'^ (B.8) 

with the decay constant 7 = max||A|,A G a{K)Y While this provides the desired 
exponential decay, the constant in front of 7" might be large. In particular, one 
would like to have a bound which is expressed directly in terms of the initial u and 
not in the expansion coefficients cx. 

If in (B.l) u is assumed to be continuous, then the dual space is the space of 
finite complex measures over [0,r]. After one application of the dual operator K* 
the measure has a continuous part u{t)dt, < t < t, and a (5-function at r. Thus 
the natural space for K^^ is C © L^([0,r],dt) = "H. For (c, -u) G 7/, c corresponds 
to the complex weight of the 5-function at r and u to u{t), < t < r. With this 
representation the operator can be written as 

where /3 = e-("+^^)^, g{t) = e-("+''^)*, f{t) = e-^"-^^)^^-*), and the kernel of A, 
A{t, s) = xi-s < t)ae-(°+''^)(*~"). 

To estimate (-ft't^)" for large n we use the resolvent of K^. As verified directly 
by integration, one finds 

{q - Ay\t, s) = q-^6{t -s) + q-\{s < t)«e-("+'")(*-^)e«"'"(*-^) (B.IO) 

for all g G C \ {0}. K^^ is a one-dimensional perturbation of A. Hence 





(q-A) 



-1 



where 

D{q) = q-P-{f\{q-A)-'g). (B.12) 
The scalar product is obtained to 

- A)-^g) = e-(°+i'^)V~'"" - e-("+^'^)^ (B.13) 
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and hence 



D{q) = q 



(B.14) 



Setting g = e^ one observes that D[q) = is equivalent to (B.4) and (B.ll) holds 
whenever g G C \ a^K^^). 

From dB.lOD, using g{t), f{t) as below (|R9|), 



+aq (1 — 2iujqa 



iq-A)-'g{t)=q-'e 



-lN-l„-(a-ia;)(T-t)„g-ia(r-t) 



m, (B.15) 

(B.16) 



The poles of the resolvent define the spectrum a{K^). There is a pole at g = 0, 
which does not contribute to the spectral representation of K^. There is an apparent 
pole at g = /3, which can be lifted since the numerator vanishes at g = /3. Finally 
there are the poles resulting from the zeros of D, 



a{K^) = {\\D{\) = 0} 



(B.17) 



which has been shown already to be in agreement with (B.4). The residue is com- 
puted from 

D'{\) = 1 + t\-K (B.18) 



Altogether, with C\ a small contour around A, one obtains 



Px = i dq{q-K^)-^ = {l + aT\-^)-^\g^){h\ 

JCx 



for A G (j{K^) with 

\gx) = {l,gx{t)), |/,) = (l,/,(t)). 
With these preparations we obtain 

K^= Yl K^ + c^r\-'r'\9x){fx\- 



(B.19) 



(B.20) 



(B.21^ 



There are no Jordan block's since the projectors in (B.21) are one-dimensional. As 
a consequence 

iK^r= A"(l + arA-i)-i|(7A)(/A|. (B.22) 

For the particular case a; = 0, 1 + ar > 0, we conclude that 
hm (KoTuit) = 9iit){fi\n) 

n— >oo 

= ( 1 + arj \u{T) + a [ dsu{s)) , (B.23) 
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as claimed in (1.10). 



To find the rate of convergence, under the same conditions, we subtract Pi and 

use 

||i^o"-^i||< E |Ari(l + arA-i)|||^A||||/A||. (B.24) 



Computing the norms in (B.24) one uses 



|A + «r|-i|A| \\gx\\ \\fx\\ = |A + aT\-'\X\ (l + |A + A* - 2\X\X\l - \X\^)) < 

(B.25) 

for all A G d-{K). The last inequality can be checked directly using that A = e^ 



satisfies (2.15). Hence 



In particular. 



\K -Pi\\< Yl l^r(«^)"' = (Amax)"(ar)-^ . (B.26) 



ds(iro"w(s)-M)'<(A„,ax)'"(ar)-M dsu{s)\ (B.27) 

^0 



which is a more explicit bound than (B.8). 
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